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The basic entity of nonlinear interaction in Navier-Stokes and the Magnetohydrodynamic (MHD) 
equations is a wavenumber triad (k,p,q) satisfying k + p + q = 0. The expression for the combined 
energy transfer from two of these wavenumbers to the third wavenumber is known. In this paper 
we introduce the idea of an effective energy transfer between a pair of modes by the mediation 
of the third mode, and find an expression for it. Then we apply this formalism to compute the 
energy transfer in the quasi-steady-state of two-dimensional MHD turbulence with large-scale kinetic 
forcing. The computation of energy fluxes and the energy transfer between different wavenumber 
shells is done using the data generated by the pseudo-spectral direct numerical simulation. The 
picture of energy flux that emerges is quite complex — there is a forward cascade of magnetic energy, 
an inverse cascade of kinetic energy, a flux of energy from the kinetic to the magnetic field, and a 
reverse flux which transfers the energy back to the kinetic from the magnetic. The energy transfer 
between different wavenumber shells is also complex — local and nonlocal transfers often possess 
opposing features, i.e., energy transfer between some wavenumber shells occurs from kinetic to 
magnetic, and between other wavenumber shells this transfer is reversed. The net transfer of energy 
is from kinetic to magnetic. The results obtained from the studies of flux and shell-to-shell energy 
transfer are consistent with each other. 



PACS: 47.65. +a,47.27-i,47.11.+j 

I. INTRODUCTION 

In fluid and MHD turbulence, eddies of various sizes interact amongst themselves; energy is exchanged among 
them in this process. These interactions arise due to the nonlinearity present in these systems. The fundamental 
interactions involve wavenumber triads (k,p,q) satisfying k + p + q = 0. The combined energy transfer computation 
to a mode from the other two modes of a triad has generally been considered to be fundamental. This formalism has 
played an important role in furthering our understanding of interactions in fluid turbulence. In this paper we present 
a new scheme to calculate the energy transfer rate between two modes in a triad interaction in magnetohydrodynamic 
(MHD) turbulence — we call it the mode-to-mode transfer. Using our scheme we calculate cascade rates and energy 
transfer rates between two wavenumber shells. 

The energy spectra and cascade rates are quantities of interest in the statistical theory of MHD turbulence. Recent 
numerical (l|-^| and theoretical work indicate that the energy spectrum of total energy is proportional to fc~ 5 / 3 , 
similar to that found in fluid turbulence. Sridhar and Goldreich j?J and Goldreich and Sridhar B studied the 

energy spectrum in presence of strong magnetic field and found it to be proportional k ± 5 ^ 3 (kj_ is the perpendicular 
component of the wavevector). Most of the work on MHD turbulence focus on the energy spectrum of either total 
energy or on the spectra of Elsasser variables (u ± b, where u and b are velocity and magnetic fields respectively). 
In this paper we will focus on various energy cascade rates of MHD turbulence. The magnetic energy (denoted by 
ME) of a mode evolves due to two nonlinear terms, [b.(u.V)b] and [b.(b.V)u], of the MHD equations; the first term 
exchanges ME between different scales, and the second term exchanges magnetic and kinetic energy (denoted by KE) 
between different scales. The KE similarly evolves due to two nonlinear terms [u.(u.V)u] and [u.(b.V)b]; the first 
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one exchanges KE between different scales and the second exchanges energy between the magnetic and the velocity 
fields. 

Pouquet et al. |J studied the energy transfer between large and small scales using three-dimensional (3D) EDQNM 
closure calculations. For nonhclical MHD they argued that the ME cascades forward, i.e., from large scales to small 
scales. In presence of helicity, the large-scale magnetic field grows, which in turn brings equipartition the small-scale 
kinetic and magnetic energy by Alfven effect. The residual helicity, the difference between kinetic and magnetic 



helicity 1 , induces growth of large-scale magnetic energy and helicity. Pouquet and Patterson |10| arrived at a similar 
conclusion based on their numerical simulations. Their simulation results also indicated that magnetic helicity, not 
necessarily the residual helicity, is important for the growth of large-scale magnetic field. Earlier, Batchelor JllJ had 
argued for transfer between kinetic and magnetic energies at small scales; his arguments were based on an analogy 
between the magnetic field equation of MHD and the vorticity equation. 

In a 2D EDQNM study, Pouquet [Q obtained eddy viscosities for MHD turbulence. Upon injection of kinetic 
and magnetic energy, she obtained a quasi stationary-state with a direct cascade of energy to small scales together 
with an inverse cascade of mean-square magnetic vector potential. In the inverse cascade regime, she argued that 
the small-scale ME acts like a negative eddy viscosity on the large-scale ME. The inverse cascade of the mean-square 
magnetic vector potential was conjectured to arise due to the destablization of the large-scale magnetic field by the 
small-scale magnetic field. She also argued that the small-scale KE has the effect of a positive eddy viscosity on the 
large-scale ME. Ishizawa and Hattori (l3) in their EDQNM calculation obtained the eddy viscosity due to each of 
the nonlinear terms in MHD equations. They found that the eddy viscosity due to (b.V)u is positive, leading to a 
transfer of energy from large-scale magnetic field to small-scale magnetic field, apparently contradicting Pouquet 's 
work p2| . However, this discrepancy may be because Pouquet considers the regime of inverse cascade of mean-square 
vector potential, while Ishizawa and Hattori p"3| ] calculation may capture the other wavenumber range. One of the 
similarities between Pouquet jlj] and Ishizawa and Hattori's jl3| calculations is that both of them give a non-local 
energy transfer from small-scale velocity field to the large-scale velocity field. In a recent work Ishizawa and Hattori 
JuJ employed the wavelet basis to investigate energy transfer in 2D MHD turbulence. 

In this paper we have computed various cascade rates and other forms of energy transfers using 2D MHD simulations. 
It is in anticipation that some of the conclusion drawn here will shed light on the nature of interactions of both 2D 
and 3D MHD because the nature of cascade of the total energy is similar in 2D and 3D . It is known that the 
ME decays in 2D MHD turbulence However, this decay occurs only after a finite time, during which the ME is 
steady. We have calculated our fluxes and shell-to-shell energy transfer rates for this "quasi steady state" . We must 
however add that the methods used here are completely generalizable to the 3D case. We have restricted ourselves to 
2D purely because of unavailability of powerful computing resources to us. 

Most of the earlier work on energy transfer in MHD turbulence (e.g., P,|lO|]) often dealt mainly with coarse-grained 
energy transfer (between large scales and small scales) . In a recent work, Frick and Sokoloff JItJ solved a shell model of 
MHD turbulence and calculated only the KE flux between the velocity modes, and the ME flux between the magnetic 
modes. In our simulations we investigate energy transfer between u-to-u, b-to-b, and u-to-b modes. These results 
provide us with more informed picture of the physics of energy transfer in MHD turbulence. 

The paper is organized as follows: In Section || we discuss the known results of energy transfer in a MHD triad, 
and introduce a new formalism called "mode-to-mode energy transfer" . Using the formalism of the mode-to-mode 



transfer, we derive in Section [II the formulas for various energy fluxes and shell-to-shell energy transfer rates. The 



simulation methodology is discussed in Section IV. In Section M we report the numerically calculated energy fluxes 
and shell-to-shell transfer rates. The results are summarized in the last section. Appendix contains derivation of the 
formulas for the mode-to-mode transfer rates. 

II. ENERGY TRANSFERS IN A MHD TRIAD 

We will first state the known result on the energy transfers in a triad. Then we will introduce a new formalism 
called "mode-to-mode" energy transfer in turbulence. 
The incompressible MHD equations in real space are 

-^ + (u-V)u = -Vp+(b- V)b + iA7 2 u (2.1) 
ot 



1 Kinetic helicity^ J u.ujd 3 x, where w = Vxu. Magnetic current helicity = J h.jcPx, where j = V x b 
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— + (u-V)b = -(b-V)u + ryV 2 b (2.2) 



at 



V-u = (2.3) 
V • b = (2.4) 



where u and b are the velocity and magnetic fields respectively, p is the total pressure, and v and r\ are the kinematic 
viscosity and magnetic diffusivity respectively. We have assumed that the mean magnetic field is zero. In Fourier 
space, the KE (E u {k) = |u(k)| 2 /2) and the ME (E b (k) = |b(k)| 2 /2) evolution equations are §JJL|1 



' j£ " (ki <.2vk?E»Qt) = ^"(k|p,q) + 7^ b (k|p,q) ( 2 - 5 ) 



k+p+q=0 k+p+q=0 



dE b (k) 



■2 V k 2 E\k)= J2 ^ 66 (k|p,q)+ E ^Wp.*!) ( 2 - 6 ) 



k+p+q=0 k+p+q=0 

The four nonlinear terms S utl (k|p,q), 5 u6 (k|p, q), S* 6 (k|p, q), and S 6 "(k|p,q) are 



5" u (k|p, q) = r"(k| P |q) + r"(k|q|p) (2.7) 

S bb (k|p, q) = ^ b (k|p|q) + Ak|q|p) (2.8) 

S ub (k\p, q) - r 6 (k|p|q) + r b (k|q|p) (2.9) 

S bu (k|p, q) = ^(k|p|q) + ^"(k|q|p) (2.10) 



where 



r"(k|p|q) = -3([k.u(q)][u(k).u(p)]) (2.11) 

^ b (k|p|q) = -3([k.u(q)][b(k).b(p)]) (2.12) 

r b (k|p|q) = 3([k.b(q)][u(k).b(p)]) (2.13) 

^"(k|p|q) = 3([k.b(q)][b(k).u(p)]) (2.14) 

where 3 stands for the imaginary part of the argument. 

The terms 5^ a (k|p, q), where a, j3 stand for u or b, are conventionally taken to represent the nonlinear energy 
transfer from modes p and q to mode k [[) 18 19[| . In the following discussion we derive a more detailed expression 



to obtain the energy transfer between any two modes within a triad; we will refer to this as the "mode-to-mode 
transfer". We emphasize that this approach is still within the framework of the triad interaction, i.e., the triad is 
still the fundamental entity of interaction of which the mode-to-mode energy transfer is a part. An expression for the 
energy transfer between two modes of a triad by the mediation of the third mode is sought here. 

Suppose in the triad (k, p, q), ^?." u (k|p|q) denotes the KE transfer from mode p to mode k with mode q as a 
mediator, and ^£. utl (k|q|p) denotes the KE transfer from mode q to mode k with p as a mediator, then by definition, 

r U (k| P |q) + r"(k|q|p) = S""(k|p, q). (2.15) 

In addition, the energy transfer from u(k) to u(p), ^?. Mtl (k|p|q), should be equal and opposite to the transfer from 
u(p) to u(k), ^ u "(p|k|q), that is, 

# bb (k|p|q)+^ bb (p|k|q) = 0. (2.16) 

Similar equations can be written for other pair of modes, and also for b-to-b and it-to-6 transfer (see Fig. [j]). Our 
objective is to obtain an expression for T]L (k|p|q) in terms of k, p, q, u(k), u(p), u(q), b(k), b(p), b(q) . 

It is shown in the Appendix that # uu (k|p|q) = ^ ,m (k|p|q) + X A ; ^ bb (k|p|q) = ^ bb (k|p|q) + r A , and ^ b "(k|p|q) = 



5# b "(k|p|q) + Aa, where ^ Q (k|p|q)s are given by Eq. (2.11-2.14). In Fig. we illustrate all these mode-to-mode 
transfers ^?. a,3 (k|p|q). The constants Aa, ^a, have special significance, and are described below. 

We take KE transfer rate as an example. In Fig. || we observe that j?""(k|p|q) + Aa gets transferred from u(p) to 
u(k), 5#™(q|k|p) + Aa gets transferred from u(k) to u(q), and $""(p|q|k) + Aa gets transferred from u(q) to u(p). 
The quantity Aa flows along u(p) — > u(k) — > u(q) — > u(p), circulating around the entire triad without changing 
the energy of any of the velocity modes. Therefore we will call it the "circulating transfer". The circulating transfer 
lost by a mode is regained, hence, this term does not affect the overall energy input/output of a mode. Therefore 
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j£" u (k|p|q) [cf. Eq. (2.11)] is the only term which effectively participates in the energy transfer. That is why it is 
termed as the "effective mode-to-mode transfer rate" from mode u(p) to mode u(k), mediated by mode u(q). There 
are two other circulating transfers: Ya for 6-to-6 transfer, and for u-to-b transfer; they are illustrated in Fig. [| A 
derivation of the effective mode-to-mode transfer is given in Appendix. For further details, refer to Ph. D. thesis of 
Dar |0| and Dar et al. 

In summary, we have obtained a formula for "effective mode-to-mode energy transfer rates" in MHD turbulence. 
The effective KE transfer rates from mode p to k is $ utl (k|p|q), the effective ME transfer rates from mode p to k is 
$ (k|p|q), and the effective "conversion" rate of KE of mode p to ME of mode k is $ "(k|p|q). The rate $ "(k|p|q) 
is responsible for the growth of magnetic energy. In all these effective energy transfers, the mode with wavenumber q 
mediates the transfer. 

In the next section we will use our "effective mode-to-mode energy transfer rate" formalism to derive formulas for 
cascade rates and shell-to-shell energy transfer. 



III. CASCADE RATES AND SHELL-TO-SHELL ENERGY TRANSFER IN MHD TURBULENCE 

In this section we will use the formalism of effective mode-to-mode transfers to define cascade rates and shell- 
to-shell energy transfer rates in MHD turbulence. Note that the formulas for cascade rate etc. derived using our 
mode-to-mode transfer is equivalent to those derive with S(k\p, q). However, our formalism has certain advantages. 
For example, some of the cascade rates and shell-to-shell transfers (defined below) which were not accessible with 
the earlier formalisms can now be calculated using our scheme. In addition, the flux and shell-to-shell transfer 
formulas using our schemes are simpler. In this paper we shall use the term ti-sphere (shell) to denote a sphere (shell) 
wavenumber-space containing velocity modes. Similarly, the 6-sphere (shell) denote the corresponding sphere (shell) 
of magnetic modes. 



A. Energy cascade rates 

The energy cascade rate (or flux) is defined as the rate of loss of energy from a sphere in the wavenumber space to 
the modes outside. There are various types of cascade rates in MHD turbulence. We have schematically shown these 
transfers in Fig. |[ The energy transfer could take take place from inside/outside w/6-sphere to inside/outside-u/6 
sphere. In terms of <# Q/3 (k|p|q), the energy cascade rate from inside of the a-sphere of radius K to outside of the 
/3-sphere of the same radius is 

n^W= E E ^ Q ( k lp|q)- (3- 1 ) 
iki>K| P i<K 

where a and fS stand for u or b. 

There is a transfer of energy from a w-sphere of radius K to the 6-sphere of the same radius. We can obtain the 
effective flux from the w-sphere to the 6-sphere using 

K<(K)= E E ^"(k|p|q). (3.2) 

IpKif |k|<K 

Similarly, the energy flux from modes outside of the u-sphere to the modes outside of the 6-sphere can be obtained by 

K>(K)= E E ^"(k|p|q). (3.3) 

|p|>K |k|>_ftT 



The total effective flux is defined as the total energy (kinetic+magnetic) lost by the if-sphere to the modes outside, 
i.e., 

n tot (K) = H«< (K) + U b < (K) + H^< (K) + H*< (K) . (3.4) 



A schematic illustration of the effective fluxes defined in Eqs. (p.l|)-( |3.3| ) is shown in Fig. |3J. In Section V^B we will 
present the values of these fluxes calculated using numerical simulations of 2D MHD turbulence. 
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B. Shell-to-Shell energy transfer rates 



In this section we derive an expression for the energy transfer rates between two shells. Consider KE exchange 
between shell m and shell n shown in Fig. |J. For the triads shown in the figure, the mode k is in shell m, the mode 
p is in shell n, and mode q could be inside or outside the shells. In terms of the mode-to-mode transfer ^. utl (k|p|q), 
the rate of energy transfer from the n th it-shell to the m th it-shell is defined as 

kem pen 

where the k-sum is over ii-shell m and the p-sum is over it-shell n with k+p+q=0. As discussed in the last section, 
the quantity T^ uu can be written as a sum of an effective transfer ^ u "(k|p|q) and a circulating transfer Xa- We know 
from the last section that the circulating transfer does not contribute to the energy change of modes. From Fig. [I| we 
can see that Xa flows from shell m to shell n and then flows back to m indirectly through the mode q. Therefore the 
effective energy transfer rate from n th it-shell to m th it-shell is 

T r n = EE^ u ( k ip|q)- ( 3 - 6 ) 

kem pen 

A general formula for shell-to-shell transfer from m th a-shell to n th /3-shell is 

kem pen 

Earlier Domaradzki and Rogallo |2^] had given the following formula for shell-to-shell energy transfer in terms of 
5 utl (k|p,q): 

^ = ^EE 5Wi ( k ip>q)- ( 3 - 8 ) 

kem pen 



This formula is incorrect due to the following reasons. The formula (3.8) has contributions from two types of wavenum- 
ber triads (see Fig. |^). In both the triads, the wavenumber k is located in the shell m. However, Type I triads have 
both wavenumbers p, q in the shell n, but Type II triads have only wavenumber p in the shell n. Clearly, to compute 
the energy transfer from shell n to shell m, we must include all the triad interactions of type I, but only the mode p 
to mode k energy transfer of type II triads. Inclusion of mode q to mode k will lead to incorrect estimation. This is 
the drawback of Domaradzki and Rogallo's formalism. The above inconsistency has been hinted at by Batchelor Jll]] 
and Domaradzki and Rogallo p2[ themselves. 

In summary, Batchelor's and Domaradzki and Rogallo's formalism does not yield correct shell-to-shell energy 
transfers; we have used an alternate formula [Eq. fl3.7|)] to compute them. We ha ve n umerically computed them using 



the numerical data of 2D MHD turbulence. These results, discussed in Section V C , provide important insights into 
the energy transfer process in MHD turbulence. 



IV. SIMULATION DETAILS 
A. Numerical Method 

In our simulations we use the Elsasser variable z^ = u ± b instead of u and b. The MHD equations (1) and (2) 
written in terms of z + and z~ are 
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Ot 



V) z ± = -Wp + v±V 2 z ± + v T V 2 z T 

+ {v ± /k%) V 4 z± + (v T /kl q ) V 4 z^ + F± (4.1) 



where v± = (v ± T))/2. In addition to the viscosity and the magnetic diffusivity of Eqs. (2.1) and (2.2), Eq. (4.1) 
includes hyperviscosity (v±/k 2 q ) to damp out the energy at very high wavenumbers. We choose i/ = t) = 5x 1CP B for 
runs on a grid of size 512 x 512. The parameter k eq is chosen to be 14 for all the runs. The terms F ± are the forcing 
functions. The corresponding forcing functions for the u and the b fields in terms of F ± will be F u = (F + + F~)/2 
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and F b = (F + — F~)/2. In our simulations, we do not force the magnetic field (i.e, F & = 0). Consequently, we get 
F + = F = F" = F. The absence of magnetic forcing is motivated by the dynamo (magnetic field amplification) 
calculations where only KE is forced. Some of the physical examples where only KE is forced are galactic dynamo 
and earth's magnetic field generation. Even though our simulations are two-dimensional, we have followed the similar 
line of approach as the dynamo simulations. 

We use the pseudo-spectral method [|3| to solve the above equations in a periodic box of size 2n x 2tt. In order 
to remove the aliasing errors arising in the pseudo-spectral method a square truncation is performed wherein all the 
modes with \k x \ > N/3 or \k y \ > N/3 (N = 512) are set equal to zero. The equations are time advanced using the 
second order Adam-Bashforth scheme for the convective term and the Crank-Nicholson scheme for the viscous terms. 
The time step At used for these runs is 5 x 10 -4 . All the quantities are non-dimensionalised using the initial total 
energy (one unit) and length scale of 2ir. Normalized cross-helicity, defined as 

2 f u.bd 3 x , d . 
- - J (4.2) 



c j{u 2 + b 2 )<p x y 

has the value approximately equal to 0.1 for the results discussed in the paper. However, we have carried out 
simulations up-to a c ~ 0.9 and the results obtained for the higher er c 's were found to be qualitatively similar to those 
discussed in this paper. 

At each time step we construct F as a divergenceless uncorrelated random function (V.F = 0). The cc-component 
of F is determined at every time step by 

F x = VMTe*t (4.3) 

where T 2 is equal to the average energy input rate per mode, and phase 4> is a uniformly distributed random variable 
between and 2n. The y-component of the forcing function is obtained by using the divergenceless condition, which 
yields 

F, = - (£) F x . ,4.4, 

We apply isotropic forcing over a wavenumber annulus 4 < k < 5. The value of T 2 is determined from the average 
rate of the total energy input, which in our simulation is chosen to be equal to 0.1. 



B. Numerical computation of fluxes 

To compute the fluxes we employ a method similar to that used by Domaradzki and Rogallo [^2| . We outline this 
method below using U!^(K) as an example. In Eq. (3.1) we substitute the expression for $ u0 (k\p\q) from Eq. pl|) 



Ki(K)= J2 E 3([k.b(q)][u(k).b(p)]) (4.5) 

|k|>if \p\<K 

A straightforward summation over k and p involves 0(N 2 ) (N = grid size) operations, and would thus involve a 
prohibitive computational cost for large N simulations. Instead, the pseudo-spectral method can be used to compute 
the above flux in O(NlogN) operations. The procedure is as follows: 
We define two "truncated" variables u > and b< as follows 



and 



b<(D)-<! b(p) if |Pl <X (A 7 ) 



Eq. (4.5) written in terms of u > and b< reads as follows 



G 



Ki(K) = E E 3 (t k - b ( k - P)] [u > (k).b<(p)]) 



(4. 



The above equation may be written as 



n fc< 



(K) 



•(k)5>(k-p)6<(p) 



(4.9) 



The p summation in the above equation can be recognized as a convolution sum. The right hand side of Eq. ( |4.9| ) 
can be efficiently evaluated by the pseudo-spectral method using the truncated variables u> and h < . This procedure 
has to be repeated for every value of K for which the flux needs to be computed. The rest of the fluxes and transfer 
rates defined in Section III are similarly computed. 

In the following section we describe the results of our simulations. 



V. SIMULATION RESULTS 



A. Generation of Steady State 



We compute the fluxes and the shell-to-shell transfers for steady state. However, the computational time required 
to obtain a statistically steady state on a grid of size 512 2 is large. So to obtain a steady state in simulations on this 
grid, we proceed in stages. First, we run on a grid of size 64 2 till a steady state is achieved. This steady state field is 
then used as the initial condition to achieve a steady-state for 128 2 , and similarly to the grid of sizes 256 2 and 512 2 . 

Theoretically, the magnetic energy in two-dimensional MHD decays in the long run even with steady KE forcing 
p6| . However, we find that the ME remains steady for sufficiently long time before it starts to decay. For 128 2 
simulation, the decay of ME starts only after t = 25 — 30 time units, and for 512 2 , the decay starts much later. We 
term this region of constant ME as "quasi steady-state" . In Fig. || we show both KE and ME for this state. Here 
the Alfven ratio (KE/ME) fluctuates between the values 0.4 and 0.56, hence ME dominates KE. The flux and the 
shell-to-shell transfer rates reported in this paper is averaged over this quasi steady-state. The averaging is done once 
in every unit of non-dimensional time over 15 time units. The averaged spectra of KE and ME over this period is 
shown in Fig. |^. For the intermediate range of wavenumbers (inertial range), the spectra show a powerlaw behaviour. 
The spectral indices are in the range of 1.5-1.7, hence the determination of MHD spectral index is quite difficult. 
However, from the flux studies, Verma et al. |l|] had concluded that Kolmogorov's 5/3 exponent is preferred over 
Kraichnan's 3/2 exponent. Miiller and Biskamp" and Biskamp and Miiller arrived at a similar conclusion based on 
their high resolution simulation. 



B. Flux studies 



We compute all the energy fluxes defined in subsection III A by averaging over the outputs of our numerical 
simulation. In Fig. [?] we show all the fluxes. The total flux n fot is positive indicating that there is a net loss of 
energy from the if-sphere to modes outside. In the wavenumber region 25 < K < 50, the total flux is approximately 
constant. This wavenumber region is the inertial range. 

The net transfer from KE to ME is a sum of the fluxes n^<(A'), U^<(K), U^>, and U^>(K). We observe from 
Fig. ^that the fluxes ^y(K max ), U^(K max ) and (K max ) are zero, understandably so because there is no mode 
outside this sphere of maximum radius K max . The flux Ii^(K max ) is found to be positive, indicating that there is 
a net transfer from kinetic energy to magnetic energy. 

We remind the reader that the w-modes within a wavenumber sphere are called the it-sphere, and the 6-modcs 
within the wavenumber sphere are called the 6-sphere. We find that the flux n^_5 (K) is positive — hence, KE is lost 
by a w-sphere to the corresponding 6-sphere. The flux (K) is also positive. It means that a u-sphere loses energy 
to the modes outside the 6-sphere. We find that the flux (K) is negative, which implies that the 6-sphere gains 
energy from modes outside the M-sphere. Thus, all these fluxes result in a transfer of KE to ME. However, H^(K) 
is negative, implying that there is some feedback of energy from modes outside the 6-sphere to modes outside the 
zt-sphere. The net transfer, however, is from KE to ME. We shall later show that the flux (K) plays a crucial role 
in driving the kinetic-to-kinetic flux (K) . 

The energy gained by the 6-spheres is found to be approximately constant for the range 20 < K < K max . This can 
be seen from the plot of n^S (K) + . The constancy of the flux implies that a 6-sphere of radius K and that of 
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radius of K + AK get the same amount of energy from the it-modes. Thus there is no net energy transfer from the 
it-modes into a 6-shell (of thickness AK) of the inertial range. We therefore conclude that the net energy transfer 
from it-modes to the 6-sphere occurs within the K < 20 sphere (the large length scales). 

We find that the flux n™> (K) is negative, that is the inertial range it-spheregain energy from the outside it-modes — 
this is consistent with the numerical simulations of Ishizawa and Hattori p3 . This behaviour of kinetic energy is 
known as an "inverse cascade" in literature and is reminiscent of the inverse cascade of KE in 2D fluid turbulence |19f] 
and of mean square vector potential in 2D MHD turbulence |jl2). Note, however that the KE in 2D MHD turbulence 
is not an inviscid invariant. Fig. |^ shows that the inverse cascade of kinetic energy exists in the wavenumber range 
K < 60. Even the modes that are being forced (4 < k < 5) gain energy from higher wavenumbers it-modes. The 
source of this energy is the positive flux II* > (K) from higher 6-modes to the higher it-modes. 

We also observe from Fig. ^ that there is a loss of magnetic energy from the 6-sphere to the 6-modes outside [see 
n|>(-K")], i.e., the ME cascade is forward. This result is consistent with the numerical results of Ishizawa and Hattori 
fl4l| . The flux H^5(K) is found to be constant in the inertial range. Using similar reasoning to that given above, we 
can again conclude that the net magnetic energy transfer to a wavenumber shell in the inertial range is zero. 

The above mentioned results are schematically illustrates in Fig. |§| for an inertial range wavenumber K = 20. The 
energy input due to forcing and the small inverse cascade [n">(.ff)] into the it-sphere from higher it-modes provides 
the energy input into the w-sphere. This energy is transferred into and outside the 6-sphere by H^(K) and H^(K) 
respectively, the latter transfer being the most significant of all transfers. The energy transferred into the 6-sphere 
from the it-sphere [n^5 (K)] and a small input from the modes outside the it-sphere [— II£>(if)] cascades down to 
the higher wavenumber 6-modes. In the high wavenumber 6-modes, the cascaded energy [11^ (if)] and the energy 
transferred from the it-sphere are partly dissipated and partly fed back to the high wavenumber it-modes [Il'>(Jf)]. 
The feedback to the kinetic energy is mostly dissipated, but a small inverse cascade takes some energy back into the 
it-sphere. This is the qualitative picture of the energy transfer in 2D MHD turbulence. 

The net transfer to each of the four corners of the Fig. || sum to zero within the statistical error (which is computed 
from the standard deviation of the sampled data). This is consistent with a quasi-steady-state picture. The results 
presented here for the quasi-steady-state in a forced turbulence remain qualitatively valid even for a decaying case — 
the direction of the various fluxes for the decaying case are identical to that for the forced simulation; but as the 
energy decays, the magnitudes of all the fluxes decrease. 

The fluxes give us information about the overall energy transfer from inside/outside it/6-sphere to inside/outside 
u/6-sphere. To obtain a more detailed account of the energy transfer, energy exchange between the wavenumber shells 
have been studied; these results are presented in the next subsection. 



C. Shell-to-Shell energy transfer-rate studies 

Significant details of energy transfers are revealed by calculating the shell-to-shell energy transfer rates defined 
by Eq. (O). We partition the k-space into shells at wavenumbers k n (n = 1, 2, 3, ...) = 1, 16, 19.02, 22.62, 2("+ 14 )/ 4 . 



The first shell extends from k\ — 1 to ki = 16, the second shell extends from &2 = 16 to k^ = 19.02, the m th shell 
exten ds f rom k m to kr m +i)- Thus, the effective shell-to-shell energy transfer rate from the n^a-shell to the m th /3-shell 



Eq. (3.7)] can be written as 



k m <k<k m +i k n <p<k„ +1 q 

The energy transfers (u — u, b — 6, and u — b) involving shells which are close in wavenumber space are called local 
transfers, as is the convention. Two shells are considered close if the wavenumber ratio of the larger shell to the smaller 
shell is less than 2 (in this study this would include the shells between n + 4 to n — 4 for the nth shell) . Transfers 
involving shells more distant than the above range are called non-local. Since two fields are involved in MHD, there 
are energy transfers between shells of same index. In the following discussions we will show that the energy transfer 
between u and b shells of same index plays a major role in MHD. 

In Figs. [I [l(| and |l| we plot the energy transfer rates T** n , T^ n and T^ n vs. m respectively for various values of 
n. It is evident from the figures that the transfer rates between shells in the inertial range are virtually independent 
of the individual values of the indices m and n, and only dependent on their differences. This means that the transfer 
rates in the inertial range are self-similar. The differences in T^" for various n are smaller than the standard deviation 
of the sampled data, indicating that the perceived self-similarity is statistically significant. 

In T^ n versus m plot of Fig. ^ we find that the transfer rates T^ n are negative for m < n and they are positive 
for m > n. Hence a &-shell gains energy from the 6-shells of smaller wavenumbers and loses energy to the fe-shells 
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of larger wavenumbers. Hence, ME cascades from the smaller wavenumbers to the higher wavenumbers (forward). 
Since T^ n is self-similar, the energy lost from a shell (n — An) to n is equal to the energy lost from n to the shell 
(n + An). Thus, the net magnetic energy transfer into any inertial range shell is zero. In addition, we find that the 
most significant energy transfer takes place among n — 1, n and n + 1 shells. Hence, b — b energy transfer is local. 

In T^ n versus to of Fig. we find that the most dominant transfers are from the n th it-shell to (to = n± l) th shells. 
From the sign of these transfers we see that KE is gained from (n — l) th shell and lost to (n + l) th shell. This means 
that the local transfers from the adjacent shells result in a forward cascade of KE towards the large wavenumbers. The 
transfers from other shells are largely negligible except for the transfer to the shell m — 1 (shown boxed on the left 
of the figure), which represents a loss of KE from high wavenumber modes to the to = 1 shell. This nonlocal tra nsfer 
(to the first shell) produces an inverse cascade of KE to the small wavenumbers, observed earlier in Section V B . 
Thus u — u energy transfer has a forward energy cascade involving local interactions, and an inverse cascade involving 
nonlocal interactions with the first it-shell. The latter dominates the former in 2D turbulence. 

We now discuss our simulation results for the energy transfer rates between it-shells and &-shells. In Fig. |ll| we find 
that T^ n is positive for all to, except for m — n — 1 and n. This implies that the n th it-shell loses energy to all but 
(n — l) th and n th 6-shells. However, the energy gained by n th and (n — l) th 6-shell is larger than the total loss of 
energy by the it-shell thr ough all other transfers. Consequently, there is a net gain of energy by the w-shells in the 
inertial range. In Section VB it was shown that outside 6-sphere loses energy to outside u-sphere; this energy flux is 
primarily due to the above shell-to-shell transfer. 

From Fig. [n] we also observe that T^ n for to < n is smaller than that for m > n. Consequently, there is a net 
transfer of energy from a it-shell to higher wavenumber &-shells. From Fig. |ll] we also find that the energy transfer rate 
T±™ from the n th it-shell to the 1 st 6-shell is positive, implying that the 1 st &-shell gains energy from the it-shells — this 
is a no nloca l transfer of energy from the large it-modes to the small 6-modes. The negative flux n„> observed in 
Section VB (see Fig. ||) is due to this nonlocal transfer. Hence, the results on flux and the transfer rates are consistent 
with each other. 

In Fig. [jj] we show the transfer rates T%? n from the n = 1 it-shell to all the &-shells. Comparing the magnitudes 
of T°£ n we notice that the energy transfer rate from the 1 st it-shell dominates the transfers from all other it-shells. 
Hence, there is a large amount of non-local transfer from the first it-shell to the larger ^-shells, reminiscent of positive 
and large n"5. In Fig. |l3| we plot the net KE transferred into a 6-shell (=Xm^mn) (° m the figure), and also KE 
transferred from all it-shells except the first (+ in the figure) . The figure clearly shows the net energy transferred into 
the inertial range 6-shell to be nearly zero, but the KE transferred from the shells n > 2 is negative. This implies 
that the inertial range 6-shells lose energy to the it-modes of shells n > 2, but they gain roughly the equal amount of 
energy from n = 1 it-shell through nonlocal interaction. The overall transfer to the inerti al ra nge 6-shells is negligible. 
This result is consistent with the picture obtained from the flux arguments (Subsection VB, Para 4). 

We schematically illustrate in Fig. |l4| the energy transfers between shells. In this figure we show directions of the 
most significant energy transfers. The arrows indicate the directions of the transfers, and thickness of the arrows 
indicates the approximate relative magnitudes. Since the local transfer rates are self-similar, the transfers from any 
other shells in the inertial range will also show the same pattern. An inertial range &-shell gains energy from the 
smaller 6-shells (T~~ n ) and the smaller it-shells (T^) through local and non-local transfers respectively. The energy 
gained by a 6-shell from the smaller 6-shells is exclusively lost to the larger 6-shells, and the energy gained from the 
it-shells is mainly lost to the same index it-shell by T^. Regarding the velocity modes, an inertial range it-shell gains 
energy from smaller it-shells by local transfers, and from the 6-shell of same shell index. An inertial range it-shell loses 
energy to higher 6-shells (locally) and to the first u-shell and the first 6-shell (nonlocally) . As illustrated in Fig. [l4|, 
there is a transfer of energy from the first it-shell to the first 6-shell. This is one of the most significant gain of magnetic 
energy from the kinetic energy. The first it-shell also loses energy to the higher 6-shell by nonlocal interactions. 

To summarize, in the inertial range the forward magnetic energy transfer and the forward kinetic energy transfers 
are local in nature. There is also energy transfer from 6-modes of shell index n to the it-modes of same indexed shell. 
There is nonlocal transfer from inertial range it-shells to the first 6-shell, and from the first it-shell to the 6-shells. 
There is also a significant energy transfer from the first it-shell to the first 6-shell. 

In this section we described various cascade rates (fluxes) and shell-to-shell energy transfer rates. It is clear that 
the complete picture is quite complex. 



VI. CONCLUSION 



In literature we find a description of energy transfer rates from two modes in a triad to the third mode, i.e., from 
u(p) and u(q) to u(k). Here we have constructed new formulae to describe energy transfer rates between a pair 
of modes in the triad (mode-to-mode transfer), say from u(p) to u(k), where the third mode in the triad acts as a 







mediator in the transfer process. 

In our formalism the mode-to-mode energy transfer can be expressed as a combination of an "effective transfer" 
and a "circulating transfer" . The circulating transfer will not result in a change of modal energy, since the amount 
of circulating transfer gained by mode k from mode p is also lost by mode k to mode q (see Fig. ||). Only the 
effective transfer is responsible for modal energy change. As the circulating transfer does not have any observable 
effect on the energy of the modes, it may be correct to ignore it from the study of energy transfer. Using the notion 
of effective modc-to-mode transfer and the circulating transfer, we defined effective cascade rates and shell-to-shell 



energy transfer rates in MHD turbulence (see formulas of Section HI). Some of these energy transfers can not be 
calculated using the earlier formalisms, but they are accessible through our "mode-to-mode" energy transfer scheme. 

Some of our formulas are directly applicable to fluid turbulence. For example, the formulas for mode-to-mode 
energy transfer ^" u (k|p|q), the flux n™>, and the shell-to-shell energy transfer T^ n can be used for fluid turbulence 
calculations without any alteration. 

We investigated the features of kinetic and magnetic energy transfer between various scales at low values of cross 
helicity in a quasi-steady state of forced 2D MHD turbulence. Several interesting observations were made in our 
simulations. A summary of the results is given below. For the following discussion refer to Figs. |^ and [w]. 

1. There is a net transfer of energy from the kinetic to the magnetic. 

2. There is an energy transfer to the large-scale magnetic field from the large-scale velocity field and also 
from the small-scale velocity field (±1^5,). The former transfer is of a greater magnitude than the latter, hence 
the magnetic field enhancement is primarily caused by a transfer from the M-sphere to the 6-sphere. Indeed the 
first few 6-modes get most of this energy. Other significant transfers between the velocity and magnetic field 
are from the large-scale velocity field to the small-scale magnetic field (n^T ) and an interesting reverse transfer 
from the small-scale magnetic field to the small-scale velocity field (II£> ) . 

3. The picture obtained from the shell-to-shell it-to-6 transfer is more complex, but is consistent with the flux 
picture discussed in item 2. The first u-shell loses a large amount of energy to the inertial range 6-shells. On 
the contrary, the first 6-shell gains energy from the inertial range u-shells. Both these processes are nonlocal. 
Also, the positivity of IK^ and negativity of is mainly due to the above two shell-to-shell energy transfers. 
In addition to these interactions, an inertial range nth u-shell loses energy to all but (n — l) th and n th shells. 
The quantity T°£ n for m < n is smaller than that for m > n, hence the inertial range w-shells mainly lose energy 
to the higher wavenumber 6-shells. The negative n^ - is primarily due to the energy transferred from the nth 
&-shell to nth w-shell. 

4. There is a forward cascade of ME towards the small scales. This is consistent with other recent numerical 
simulations EDQNM closure calculations also yield a ME transfer to the small-scales 

5. The shell-to-shell transfer between b-to-b modes shows that the most significant interaction of n-th b shell takes 
place with (n ± l)th b shells; the (n — l)th 6-shell loses energy to the nth 6-shell, which in turn loses energy to 
the (n -I- l)th shell. Hence, the b-to-b energy transfer is forward and local. 

6. There is an inverse cascade in the velocity field — this is consistent with the observation of Ishizawa and Hattori 
arising from numerical simulations |L4| and is also consistent with the EDQNM closure calculations fll2| , |l3| | . 
This inverse cascade of KE is driven by the reverse transfer of energy from magnetic to the velocity field at the 
small scales. Although the flux study points to an inverse cascade of kinetic energy, the shell-to-shell energy 
transfer rates reveal the following feature. There exists both an inverse and a forward transfer of kinetic energy. 
The inverse transfer is primarily nonlocal, coming from the large wavenumber u-modes to the small wavenumber 
u- modes. The forward cascade is local, i.e., between same sized eddies in the inertial range. 

7. The shell-to-shell energy transfers between the inertial range shells statistically depend only on the difference 
of the shell indices. Hence, the inertial range interactions can be said to be self similar. 

As pointed out in the first section, Pouquet Jjjj observed (based on her 2D EDQNM studyj a forward cascade 
of energy and an inverse cascade of mean-square vector potential, and Ishizawa and Hattori |l3j proposed that the 
large-scale magnetic energy is enhanced due to energy transfer from the small-scale velocity field to the large-scale 
magnetic field. We also find these transfers in our 2D simulation, but the dominant transfer to the large-scale ME is 
from the large-scale KE. Our simulation can not probe the inverse cascade region referred to by Pouquet jl2| because 
we force the small wavenumber modes, so there is no scope of inverse cascade in our simulation. 

Pouquet and Patterson proposed in context of 3D turbulence that there an inverse cascade of ME from small 
length-scales to large length-scales in presence of helicity. Since there is no helicity present in 2D MHD, helical MHD 
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is inaccessible to our simulation. Still we observe some very striking features in our simulation. For example, our 
numerical result regarding forward cascade of ME (11^5 > 0) is in agreement with the forward cascade of ME obtained 
by the perturbative calculation of nonhelical MHD Similar analytic calculation in presence of helicity is in 

progress. Full three-dimensional simulation and analytic work in this direction will be useful in resolving various 
issues in the magnetic field generation problem. 

The fluxes and the transfer rates discussed here and in the analytic paper by Verma dJM] com d find applications 
in the dynamo problem. In astrophysical objects, like galaxies and the sun, the magnetic field is thought to have 
arisen due to the amplification of a seed magnetic field. Our study has been performed over a quasi-steady state 
with a low Alfven ratio. In order to understand the build-up of magnetic energy starting from a seed value it is 
important to perform a similar study at high Alfven ratio. In some of the popular models, like the a-dynamo ||l6| , the 
mean magnetic field gets amplified in presence of helical fluctuations. It must be borne in mind that our calculations 
are two-dimensional and devoid of magnetic helicity and kinetic helicity A three-dimensional calculation (with the 
inclusion of magnetic and kinetic hclicities) of various fluxes and shell-to-shell energy transfer rates will yield important 
insights into some unresolved issues concerning enhancement of magnetic energy. 

The fluxes also find important applications in various phenomenological studies. For example, Verma et al. psj j 
estimated the turbulent dissipation rates in the solar wind and obtained the temperature variation of the solar wind 
as a function of solar distance. The various cascade rates discussed here could be useful for various astrophysical 
studies. For example, 11^5 and IK^ can be used for studying the variation of ta of the solar wind. 

The physics of MHD turbulence is still unclear. The studies of various fluxes and transfer rates shed light at various 
aspects which will help us in getting a better understanding of MHD turbulence. 
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APPENDIX A: DERIVATION OF MODE-TO-MODE ENERGY TRANSFERS IN MHD TURBULENCE 

In this appendix we will derive the result stated in Section || that the "effective " energy transfer rate from mode p 
to mode k ("^"^(klplq)) is equal to ^ Q,3 (k|p|q). We will find in detail the mode-to-mode energy transfer from u(p) 
to u(k) within the triad (k,p,q). The other mode-to-mode energy transfer from b(p) to b(k) and from u(p) to b(k) 
can be obtained by following the same procedure. In all these mode-to-mode transfers, mode q acts as a mediator. 



1. Definition of Mode-to-Mode transfer in a triad 



In the beginning we limit ourselves to u to u transfer in a triad (k, p, q). Let the quantity 1/t (k|p|q) denote 
the energy transferred from mode p to mode k with mode q playing the role of a mediator (see Fig. ny. We wish to 
obtain an expression for 1^ uu . 

The Tjt uu should satisfy the following relationships : 

1. The sum of ^. Mtl (k|p|q) and 7?" M (k|q|p), which represent energy transfer from mode p to mode k and from 
mode q to mode k respec tive ly, should be equal to the total energy transferred to mode k from modes p and q, 



i.e., 5""(k|p,q) [see Eq. Thus, 

^(klplq) + r"(k|q|p) = S«"(k|p, q), (Al) 

Similarly, 

^(plklq) + ^"(Plqlk) = S""(p|k, q), (A2) 
^(qlklp) + ^"(qlPlk) = ^"(q|k, p). (A3) 
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2. The energy transferred from mode p to mode k, i.e., ^2.""(k|p|q), will be equal and opposite to the energy 
transferred from mode k to mode p i.e., ^. utl (p|k|q). Thus, 

r"(k|p|q)+r tl (p|k|q)=0, (A4) 

r"(k|q|p)+r tl (q|k| P )=0, (A5) 

r"(p|q|k)+r U (q|p|k)=0. (A6) 

These are six equations with six unknowns. However, the value of the determinant formed from the Eqs. ( |Al|) -([A6|) 
is zero. Therefore we cannot find unique 1jt wu s given just these equations. 



2. Solutions of equations of mode-to-mode transfer 



Consider 



r U (k|p|q)^-3([k.u(q)][u(k).u(p)]) 



(A7) 



By definitions given in Eqs. (2.7), ^""s satisfies Eqs. ( Al-|A3|). Using the triad relations hip k + p + q = 0, and the 
incompressibility constraint [k.u(k) =0], it can be seen that $ uu s, etc. also satisfy Eqs. (A4-A6). Hence, the set of 
$ uw s are one instance of the ^. lm s, i.e., ffi""( k|p|q) = <$ mi (k|p|q). However, this is not a unique solution because the 
determinant formed from the Eqs. ( Al)-([A6|) is zero. 

If another solution ^. utl (k|p|q) differs from ^""(k|p|q) by an arbitrary function Xa, i.e., ^""(klplq) = <# Mtl (k|p|q) + 
Xa, then by inspection we can easily see that Xa gets added to other two mode-to-mode transfer as well. This 
conclusion is depicted in Fig. ^. As shown in the figure, Xa simply "circulates" between the modes, hence it does 
not appear in the effective energy transfer. Note that Xa could depend upon the wavenumber triad (k, p, q,) and 
the Fourier components u(k), u(p), u(q). It may be possible to determine Xa by imposing additional constraints 
of rotational invariance, Galilean invariance, finiteness, etc. [pof . We discuss in Section [n] that the circulating transfer 
plays no role in modal energy change, hence we will not dwell upon the problem of its determination. In Section III 



we also show that the cascade rates and the shell-to-shell transfers can be written in terms of effective mode-to-mode 
transfer rates ($ a ). 

We can similarly show that (1) the ME transfer rate % from the mode b(p) to b(k) mediated by mode u(q) is 
given by 



% hh = $ 00 + Y L 



bb 



(A8) 



and the conversion rate of KE from the mode u(p) to ME of mode b(k) mediated by mode b(q), Ijl " is given by 



b u 



Za- 



(A9) 



The terms lfi lh and 1^ bu are given by Eqs. ( 2.12j ) and (2.14) respectively. The quantities Ya and Za are the "circulating 
transfers". As shown in Fig. ||, Ya circulates around b(k) — > b(p) — > b(q) — > b(k), while Za circulates around 
u(p) — > b(k) — > u(q) — > b(p) — > u(k) — * b(q) — ► u(p). The details of the above derivation can be found in pO|,pH. 
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FIGURE CAPTIONS 



Fig. 1: The mode-to-mode energy transfers sought to be determined. 



Fig. 2: The energy transfer between a pair of modes in a triad can be expressed as a sum of a circulating transfer and 
an effective mode-to-modc transfer. 



Fig. 3: Various energy fluxes defined in this section. The arrows indicate the directions of energy transfers cor- 
responding to a positive flux. These fluxes are computed using data of numerical simulations and are shown in 
Fig. 8. 



Fig. 4: The circulating transfer and effective mode-to-mode transfer from shell n to shell m with the arrows pointing 
towards the direction of energy transfer. 



Fig. 5: Evolution of the KE and the ME for the simulations on grid sizes 512 2 and 128 2 . For 512 2 a quasi-steady 
magnetic energy is obtained over the period of the simulation. It is demonstrated in the 128 2 simulation that a 
quasi-steady magnetic energy eventually decays. 



Fig. 6: The averaged KE spectrum E u (k) and ME spectrum E b (k): Averaging is done over 15 time units. 
Fig. 7: The plots of various fluxes vs. wavenumbers. 



Fig. 8: The schematic illustration of numerically calculated values of fluxes averaged over 15 time units. The values 
shown here are for K = 20, a wavenumber within the inertial range. 



Fig. 9: The energy transfer rate T^ n from the n th fe-shell to the m th 6-shell. The loss of energy from the n th 6-shell 
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to the m th 6-shell is defined to be positive. 



Fig. 10: The energy transfer rate T^™ from the n th M-shell to the m th M-shell. The boxed points represent energy 
transfer from the n th it-shell to the 1 st M-shell. 



Fig. 11: The energy transfer rate T%? n from the n th M-shell to the m th 6-shell. 



Fig. 12: The energy transfer rate from the I s * M-shell to the 6-shells. 



Fig. 13: The diamonds (o) represent the net energy transfer into a 6-shell from all the M-shells. The pluses (+) 
represent the energy transfer into a 6-shell from shells n>2. 



Fig. 14: A schematic representation of the direction and the magnitude of energy transfer between M-shells and b- 
shells. The relative magnitudes of the different transfers have been represented by the thickness of the arrows. The 
non-local transfers with the I s * shell have been shown by dashed lines. 
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